Section: Software

FVforBlochTorrey

Participant : Jing-Rebecca Li [correspondant] .

Finite volume code in Fortran 90 to solve the multiple compartment Bloch Torrey equation in 2D and 3D to simulate the bulk magnetization of a sample under the influence of a diffusion gradient. We couple a mass-conserving finite volume discretization in space with a stable time discretization using an explicit Runge-Kutta-Chebyshev method and we are able to solve the Bloch-Torrey PDE in multiple compartments for an arbitrary diffusion sequence with reasonable accuracy for moderately complicated geometries in computational time that is on the order of tens of minutes per bvalue on a laptop computer.

See also the web page http://www.cmap.polytechnique.fr/~jingrebeccali/ .